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ABSTRACT 

We undertake a quantitative investigation, using Monte Carlo simulations, of the 
amount by which quasars are expected to exceed radio galaxies in optical luminosity 
in the context of the 'receding torus' model. We compare these simulations with the 
known behaviour of the [O III] A5007 and [O II] A3727 emission lines and conclude that 
[O III] is the better indicator of the strength of the underlying non-stellar continuum. 
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1 INTRODUCTION 

It is widely believed that radio-loud quasars and radio galax- 
ies differ from each other only in terms of the angle be- 
tween their radio axes and the line of sight (Scheuer 1987; 
Barthel 1989). Quasars are observed fairly close to the line 
of sight (< 45°) and therefore frequently exhibit the effects 
of beaming, such as superluminal motion and luminous flat- 
spectrum radio cores, whereas radio galaxies are observed 
with their axes close to the plane of the sky, and do not show 
these effects. However, quasars also possess a luminous non- 
stellar optical continuum and broad emission lines, which 
are absent in radio galaxies. It has therefore been proposed 
that there is material around the nucleus which lies pref- 
erentially in the plane perpendicular to the radio axis and 
obscures the central regions from view in radio galaxies. Be- 
cause of the geometry of this material, it is often referred 
to as the "torus", but other geometries (e.g. a warped disk) 
are possible. Although the broad lines are hidden from di- 
rect view by this material, the narrow line region (NLR) is 
much larger in size and should be less strongly affected. It 
is therefore possible that narrow line luminosity could be an 
isotropic measure of the strength of the non-stellar contin- 
uum that is otherwise invisible in radio galaxies. 

Jackson & Browne (1990) tested this idea by comparing 
the [O III] A5007 luminosities of quasars and radio galax- 
ies with similar redshifts and extended radio luminosities. 
They reported that quasars are 5-10 times more luminous 
in this line than radio galaxies, and attributed the difference 
to higher extinction in radio galaxies. Hes, Barthel & Fos- 
bury (1993) performed a similar test using the [O II] A3727 
doublet and found that there was no measurable difference 
between the line luminosities of the two classes. Since [O II] 
is at a shorter wavelength than [O III] and would be more 
greatly affected by a foreground screen of reddening, they 



suggested that most of the [O III] emission is produced close 
to the nucleus and is still obscured by the torus in radio 
galaxies, whereas [O II] is produced farther out and is unaf- 
fected. 

Unfortunately, Jackson & Browne's analysis was bi- 
ased because, although the radio galaxies were selected from 
the 3CR-based complete sample of Laing, Riley & Longair 
(1983; also called 3CRR), the quasars were drawn from a 
number of incomplete surveys and were therefore subject 
to uncertain selection effects, most notably the tendency to 
preferentially include optically bright objects. Since line and 
continuum luminosities are very well correlated in quasars, 
their result would be biased towards finding systematically 
higher line luminosities in the quasars. In addition, Jackson 
& Browne's radio galaxy sample included some objects with 
very weak emission lines (Class B optical spectra; Hine & 
Longair 1979) that it is now believed may not belong to the 
unified scheme. An unbiased analysis significantly reduces 
the magnitude of Jackson & Browne's result, but the quasars 
are still about twice as luminous in [O III] than the radio 
galaxies. This result is also seen in the higher [O IIl]/[0 II] 
ratios in broad-line objects than in narrow-line objects (e.g. 
Saunders et al. 1989; Tadhunter et al. 1993). At higher red- 
shift, however, the [O III] luminosities of the two classes of 
object are comparable (Jackson & Rawlings 1997). 

Rawlings & Saunders (1991) note that the tendency for 
quasars to have a higher emission line luminosity than ra- 
dio galaxies could be the result of a classification bias. This 
could explain the discrepancy in [O III] luminosities between 
the two classes, but appears to run counter to the similar- 
ity in the [O II] line emission. In this Letter, we apply the 
simple receding torus model (Lawrence 1991; Hill, Goodrich 
& Depoy 1996) to provide a quantitative explanation of the 
difference in [O III] luminosities between quasars and radio 
galaxies (and its redshift dependence), and explain why a 
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Figure 1. Schematic representation of the receding torus model. 
The opening angle, 6, for the low luminosity object (top) is fairly 
small. In the high luminosity object (bottom), the inner radius 
of the torus, r, is larger, producing an increased opening angle 
provided the half-height of the torus, h, remains constant. 

similar effect is not seen in [0 II]. The free parameters in- 
volved in our explanation are constrained by observed quan- 
tities independent of the emission line luminosities. We con- 
clude that the luminosity of the [O III] A5007 line is an excel- 
lent indicator of the strength of the underlying non-stellar 
continuum. 



2 METHOD 

2.1 The receding torus model 

The receding torus model, as proposed by Lawrence (1991) 
and illustrated in Figure bL is fairly simple. Dust evaporates 
at a temperature of about 1500 K, and cannot survive closer 
to the nucleus than the radius where the temperature of the 
nuclear radiation field is hotter than this. All quasars display 
a near-infrared bump longward of 1 /jm that is believed to 
be caused by thermal emission from dust, and the constancy 
of wavelength at which this bump appears indicates that 
dust is always present at the hottest possible temperatures. 
The inner radius of the torus is therefore determined by the 
radius at which dust evaporates, which scales as r cx L 0,5 . If 
the half-height, h, of the torus is independent of the source 
luminosity, then the half opening angle = tan -1 r/h will 
be larger in the more luminous objects. Hill et al. (1996) 
find that this model fits the observed nuclear extinctions in 
radio galaxies very well. 

We are more likely to be viewing the luminous objects 
within their cones, and are consequently more likely to clas- 
sify them as quasars. Qualitatively then, this simple model 
could explain the tendency for quasars to have a higher 
[O III] luminosity than radio galaxies. We therefore under- 
take a quantitative investigation using Monte Carlo simula- 
tions. 

2.2 Simulations 

We consider an ensemble of objects with the same redshift 
and radio power, which could reasonably be assumed to have 
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Figure 2. Observed dispersion in a radio-selected quasar sample 
as a function of the dispersion in the parent AGN population. 8o 
is the half opening angle for the mean luminosity, Lq, and the 
solid lines are for 6 = 30°, 40°, 50°, 60°, 70°, 80°. The dashed 
line is the equality relationship. 

similar torus properties (e.g. torus height). Their intrin- 
sic optical-ultraviolet luminosities will have a well-defined 
mean, Lo, determined by the strong correlation between op- 
tical and radio luminosity, but will show scatter due to the 
intrinsic dispersion in this relationship. We assume this dis- 
persion to be Gaussian in log L, with standard deviation a. 
We define #o to be the half opening angle for an object with 
the mean luminosity, Lq. An object with luminosity L with 
therefore have an half opening angle, 

6> = tan" 1 [(L/L )° ,5 tan6io], 

and the probability that such an object will be observed as 
a radio galaxy (i.e. our line of sight is outside the opening 
angle of the cone) is 

p RG = cos 6 = (1 + tan 2 6»)~ - 5 . 

Obviously, jjqso = 1- JJrg • 

For our ensemble of objects, we can therefore determine 
the ratio of the mean optical-UV luminosities of quasars 
and radio galaxies (we evaluate the mean in log space), for 
a given a and Qq. However, there exist observational data 
which allow us to fix these free parameters. First, Serjeant 
et al. (1998) measure the scatter about the mean optical- 
radio luminosity correlation for steep-spectrum quasars to 
be 1.5 magnitudes (0.6 dex). This will be slightly smaller 
than a, since there should exist radio galaxies whose in- 
trinsic optical luminosities are on average fainter than those 
of the quasars. However, Serjeant et al.'s sample was dom- 
inated by distant (z > 1) quasars whose high luminosities 
would imply a large opening angle for the torus. At large 
opening angles, most objects are observed as quasars, and 
therefore the dispersion in optical magnitude in a quasar 
sample is close to the intrinsic dispersion of the parent AGN 
population (Figure ^). We thus adopt a = 0.6. The second 
constraint is the observed quasar (or radio galaxy) fraction; 
we must constrain 6o such that the fraction in our simulation 
matches that in a radio-selected sample. 

Following the reasoning of Laing et al. (1994) and oth- 
ers, when comparing our results to observed quantities, we 
consider only those radio galaxies with prominent emission 



© 0000 RAS, MNRAS 000, 000-000 



A new look at the isotropy of narrow line emission in extragalactic radio sources 3 




30 40 50 60 70 



Figure 3. Observed radio galaxy fraction as a function of the 
torus half opening angle, 9q, for an object with the mean lumi- 
nosity, Lq. The lines are for a = 0.3, 0.6, 0.9, 1.2. 

lines. At 2 < 0.43 and h < a < 13 h , where there exist high- 
quality optical spectra allowing classification on the basis 
of broad Ha (Laing et al. 1994), there are 15 such radio 
galaxies and 9 quasars. The observed radio galaxy fraction 
is therefore 0.625, corresponding to 80 = 50° for a — 0.6 
(Figure ^; the value is extremely insensitive to the choice of 
a). 



3 RESULTS 

3.1 Relative quasar/radio galaxy continuum and 
emission line luminosities 

Figure ^ illustrates the amount by which quasars are ex- 
pected to be overluminous compared to radio galaxies as a 
function of 9q and <r, and for the values determined previ- 
ously, we find (Lqso) / (Lug) = 1-8. This is effectively the 
same as the difference in [O III] luminosities for the two 
groups of objects, which could therefore be explained if the 
[O III] emission line luminosity is directly proportional to 
the ionizing luminosity. Its constancy of equivalent width 
over approximately three orders of magnitude in luminosity 
(Miller et al. 1992) indicates that it is. 

Why then is a similar effect not seen in the [O II] lumi- 
nosity, since both emission lines correlate well with optical 
luminosity over several orders of magnitude in quasars? We 
believe that the answer is once again due to a bias intro- 
duced by scatter. Although the effective ionization parame- 
ter of the NLR does not correlate with luminosity (the range 
of quasar luminosities is far larger than the observed range 
of U, -3 < log (7 < -1; e.g. Saunders et al. 1989), the 
most luminous objects at a given redshift will also have the 
highest ionization parameters. Figure ^ shows the depen- 
dence of [O III] and [O II] luminosity as a function of ion- 
ization parameter for an optically-thick plane-parallel slab. 
This was calculated with the photoionization code CLOUDY 
Version 84.12a (Ferland 1993) using the AGN continuum de- 
scribed in Mathews & Ferland (1987), although effectively 
identical results are obtained for power laws with a range 
of spectral index. Over the range of effective ionization pa- 
rameter typical of the NLR, the [O II] luminosity is only 
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Figure 4. Ratio of mean quasar luminosity to mean radio galaxy 
luminosity as a function of 8q for various value of a, the dispersion 
in the optical-radio luminosity correlation of the parent popula- 
tion. The lines run from a = 0.3 (bottom) to cr = 1.2 (top) in 
steps of 0.1. The bold line is for <r = 0.6. 
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Figure 5. Luminosity of the [O III] (solid lines) and [O II] 
(dashed lines) lines as a function of ionization paramater, U for an 
optically-thick plane-parallel slab. Each set of curves is for values 
of the hydrogen number density njj = 10 2 , 10 3 , 10 4 cm - 3 . 

very weakly dependent on the strength of the ionizing con- 
tinuum. Whereas a hundredfold increase in U increases the 
luminosity of [O III] by a factor of ~ 200, the [O II] lu- 
minosity becomes only ~ 4 times larger. Although a single 
plane-parallel slab is obviously an overly simplistic model 
for the NLR, there is very little density dependence in the 
way the line luminosities vary with ionization parameter. 

3.2 Redshift dependence 

The lessening of the effect at higher redshift is easy to under- 
stand, since in the flux-limited 3CR sample this corresponds 
to higher luminosity and therefore a larger mean opening 
angle. At such large angles, variations in luminosity of an 
order of magnitude or less have little effect on the angle, 
and hence the bias towards more luminous objects being 
seen as quasars is smaller. This is apparent in Figure [|. Al- 
though radio galaxies outnumber quasars 2:1 at low z, at 
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z > 1 the ratio is nearly reversed. In addition, there may 
be many objects like 3C 22 (Economou et al. 1995; Rawl- 
ings et al. 1995) or 3C 41 (Simpson, Rawlings & Lacy 1998; 
Economou et al. 1998) that possess broad Ha which is dif- 
ficult to observe from the ground at z > 1.1 due to it being 
redshifted into regions of poor atmospheric transmission or 
the H window, whose strong OH airglow lines hamper the 
detection of broad features. The observed radio galaxy frac- 
tion of ~ 40% at z > 1 should therefore be considered an 
upper limit, and it provides a lower limit to 6o •> 70°. This 
corresponds to (I/qso) / {Lb.g) < 1-7, and the ratio would 
drop below 1.5 if only one third of those objects currently 
classified as radio galaxies are really quasars. It should also 
be noted that measurement errors tend to be rather larger in 
near-infrared than in optical spectroscopy, and so the study 
of Jackson & Rawlings (1997) is less sensitive to differences 
in the two classes than are low redshift studies. 

Jackson & Rawlings (1997) do not offer any explana- 
tion as to why they fail to observe a similar effect to that of 
Jackson & Browne (1990), yet it is naturally explained in our 
model as a luminosity, rather than a redshift, dependence. 
Although our model might appear to predict that high red- 
shift samples whose luminosities are similar to those of low 
redshift 3CRR objects (e.g. a sample from the 7C catalogue) 
should also display a separation in [O III] luminosity with 
broad/narrow line class, such a prediction also depends on 
the assumption that the average torus height does not vary 
with redshift since this also controls the mean opening angle. 
We can predict that the difference between the [O III] lumi- 
nosities of radio galaxies and quasars should be a function 
of quasar fraction (Figures^ and ^), and should effectively 
vanish for samples with a large (> 80%) quasar fraction. 

3.3 Polarization properties 

Since our model requires that effectively all the [O III] emis- 
sion is seen directly, it appears to be in contradiction with 
the recent detection of polarized [O III] in radio galaxies 
by di Serego Alighieri et al. (1997). These authors calcu- 
late that 80-90% of the total [O III] emission is obscured 
in radio galaxies, but their calculations assume that the 
observed featureless continuum is entirely scattered. It is 
known that Seyfert galaxies possess an extended featureless 
continuum (FC2; Tran 1995) and it is reasonable to imag- 
ine that radio galaxies do also, thus lowering the estimate 
of the obscured fraction. In addition, the [O III] emission 
may be polarized, even if none is obscured, if the scatter- 
ing region is distributed anisotropically. Although di Serego 
Alighieri et al. disfavour this idea, it is not in conflict with 
observations, and would be expected if the scattering mate- 
rial is associated with the outflow. The fact that polarized 
[O III] emission is detected in 3C 227 despite our relatively 
clear view of the nucleus seems to indicate that significant 
obscuration is not required for the line emission to be polar- 
ized. We also note that the size of the [O IIl]-emitting region 
that must be obscured is much larger than estimates for the 
height of the torus derived from the dust evaporation radius 
and cone opening angle. While the polarization data are in- 
triguing, we therefore feel that our model is not in conflict 
with them. 



4 SUMMARY 

We have used a simple receding torus model, whose free 
parameters are constrained by observation, to show that low 
redshift 3CRR quasars should be, on average, about twice as 
luminous in their ionizing continua as radio galaxies of the 
same radio luminosity. This difference should also be seen in 
their [O III], but not their [O II], emission line luminosities, 
in agreement with observation. For samples with a higher 
quasar fraction, such as the high redshift 3CRR objects, the 
difference in ionizing luminosities between quasars and radio 
galaxies should be smaller, and there should therefore be less 
of a difference in their [O III] luminosities, again in line with 
observation. 

This model leads to the conclusion that the [O III] A5007 
emission line, and not the [O II] A3727 doublet, is an unbi- 
ased indicator of the intrinsic optical-ultraviolet luminosity 
of both quasars and radio galaxies. 
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